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Abstract: We compute rigorously the ground and equilibrium states for 
Kitaev's model in 2D, both the finite and infinite version, using an analogy 
with the ID Ising ferromagnet. Next, we investigate the structure of the 
\ reduced dynamics in the presence of thermal baths in the Markovian regime. 

Special attention is paid to the dynamics of the topological freedoms which 
■ have been proposed for storing quantum information. 

O 

1 Introduction 

Despite the enormous activity in the field of quantum information in the last 
decade, the fundamental problem of how to protect quantum information 
Q-f for macroscopic periods of time remains open. In contrast to the popular 

belief the Fault Tolerant Quantum Computation (FTQC) ideas [1, 2, 3, 4, 5] 



do not provide the final solution to this problem. The mathematical mod- 
els behind these schemes are highly oversimplified, phenomenological, and 
do not take into account the restrictions imposed by the laws of thermo- 
dynamics. As FTQC involves time-dependent control by external fields the 
rigorous analysis of the corresponding models which can be derived from first 
principles is extremely difficult. In particular, the standard approximations 
based on the rigorous versions of Markovian limits (weak coupling or low 
density [6, 7, 8]), which yield descriptions consistent with thermodynamics, 
cannot be applied, see [9, 10, 11, 12, 13] in this context. Therefore, the 
approach to quantum memories based on the idea of self-correcting systems 
is more promising for rigorous studies. It deals with iV-qubit systems with 
specially designed Hamiltonians which are supposed to stabilise an algebra of 
observables corresponding to one or few qubits, at least below some critical 
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temperature. The construction of self- correcting Hamiltonians is based on 
the properties of error correcting codes and the most advanced examples are 
2, 3 and 4D Kitaev models [14, 15] associated with toric codes. 

It is instructive to discuss first the mechanism of protection of classical 
information using phase transitions in quantum systems. In the general ap- 
proach based on the algebraic formalism for infinite systems, i.e. systems in 
the thermodynamic limit, we characterise equilibrium states by the KMS con- 
dition [16]. States that satisfy the KMS condition are known to have good 
stability properties such as passivity [17]. Typically, one expects that for 
sufficiently high temperatures there exists a unique KMS state, while below 
a certain critical temperature T c a simplex of KMS states can appear. The 
extreme states which span this simplex correspond to pure thermodynamical 
phases labelled by the values of the order parameters. They can be used to 
encode classical information. 

In a complementary dynamical picture we consider an A-particle system 
weakly interacting with a heat bath at temperature T. The finite system pos- 
sesses a unique equilibrium state given by a Gibbs density matrix. However, 
for systems exhibiting phase transitions and for T < T c the unique equi- 
librium state can be decomposed into many long-living metastable states 
corresponding to local minima of the free energy. These minima are sep- 
arated by free energy barriers of heights proportional to N. This makes 
transitions caused by thermal fluctuations highly improbable and leads to 
life-times exponentially increasing with N. Obviously, the metastable states 
should become equilibrium (KMS) states in the limit N — > oo. Although 
the equivalence of equilibrium and dynamical pictures is generally accepted 
there exist only few rigorous results that support this, e.g. the mean-field 
Curie- Weiss model [18] and Bose-Einstein condensation [19]. 

From the discussion of above it follows that extending these ideas to the 
field of quantum information in order to construct quantum memories is 
not possible. The encoding of qubits into different thermodynamical phases 
of an infinite system is impossible because of the simplex structure of the 
KMS states. A simplex is a purely classical state space and cannot even 
support a single qubit. 

Actually, in the literature one attempts to construct self-correcting Hamil- 
tonians for N qubits in such a way that their degenerated ground states 
correspond to codewords encoding states of one or few qubits [14, 15]. How- 
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ever, good codewords are locally indistinguishable, which implies that in the 
thermodynamic limit all of them merge into a single ground state of the in- 
finite system. Therefore one doesn't expect that such systems exhibit phase 
transitions at any temperature and so cannot even support a classical mem- 
ory. Such a view is consistent with Pirogov-Sinai [20] type results where one 
shows that well-separated ground states of infinite classical systems extend 
to separated phases at sufficiently low temperatures. 

A possible objection to the previous reasoning, which was presented 
in [25], is that the mathematical framework of infinite systems, used to de- 
scribe the thermodynamic limit, is too simplistic. One cannot exclude that 
observables which distinguish between the degenerated ground states of the 
finite system but disappear in the thermodynamic limit are long-lived in the 
presence of a generic heat bath at sufficiently low temperatures. This would 
imply that the expected equivalence of the infinite system equilibrium picture 
and the dynamical picture is not universal and that a new type of "purely 
dynamical phase transitions" exists. In order to prove or disprove this claim 
one should study finite supposedly self-correcting systems weakly coupled to 
generic models of heat baths and search for possibly slowly relaxing degrees 
of freedom. These could then support states of encoded qubits for times 
exponentially increasing with N. 

It is interesting to compare the problem of stable quantum memories to 
another long-standing problem of quantum mechanics: the absence of certain 
superpositions of quantum states for sufficiently large systems. Well-known 
examples are Schrodinger cat states, the appearance of a classical world, the 
problem of molecular structure, chiral molecules, deformed nuclei, . . . [21] It 
seems to be generally accepted that the main ingredient for the solution of 
this problem is the stability of quantum states of the observed system with 
respect to the coupling to an environment, sometimes identified with a mea- 
suring apparatus or even a conscious observer. As in the case of quantum 
memories two approaches have been developed, a static and a dynamical one. 
The first relies on dressing the states of the bare system. This can produce 
disjoint stable states of the dressed system when using environments with an 
infinite number of degrees of freedom and leads to unobservable superposi- 
tions, at least for coupling constant larger then a certain critical value. The 
mathematical mechanism is completely similar to that which is responsible 
for the appearance of disjoint KMS states in the thermodynamic limit. The 
dynamical approach, on the other hand, refers to decoherence phenomena 
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which transform superpositions of states into mixtures with respect to a cer- 
tain particularly stable pointer basis. The relations between both pictures 
and their ranges of applicability are still a topic of debate. Unfortunately, 
except for the spin-boson model [22] no rigorous result has been obtained. 

The purpose of this paper is to start a systematic and rigorous compari- 
son of the infinite system equilibrium picture and the finite system dynamical 
picture for models designed to support quantum memories. We restrict our- 
selves to the ID (quantum) Ising model and the 2D Kitaev model. Although 
there exist heuristic arguments that, in spite of earlier claims concerning 
the 2D case [14, 23, 24], both models cannot provide stable quantum memo- 
ries [25] it is instructive to develop the basic mathematical tools for analysing 
these relatively simple cases, hoping that they will be useful for 3D or 4D 
models. Only the last case is presently considered by some experts to be 
a proper candidate for a quantum memory [15]. One should mention that 
formally the ID quantum Ising model can be considered to be the ID Kitaev 
model and hence the parallel analysis of both could be helpful for further 
studies. 

The main simplifying feature of the Kitaev models is the absence of wave 
propagation. This leads to a discrete spectrum of the Hamiltonian with 
level spacings independent of the size of the system. This property allows 
to apply the Markovian approximation in terms of Davies weak coupling 
limit which reproduces fully the expected phenomenology of an open system 
weakly coupled to a heat bath. The Davies semigroup generators possess a 
simple local structure and behave well with respect to the thermodynamic 
limit. Moreover, for Kitaev's model we can determine the structure of the 
algebras which can encode qubits and we obtain the dynamics of the qubit 
observables. 

2 Ground and equilibrium states of ID Ising 
and 2D Kitaev models 

The physical substrate which is to carry the qubits has basic spin observ- 
ables that are attached to the sites of a chain, a ring, a lattice . . . We will 
describe their basic observables by Pauli matrices a? where j points at the 
site at which the spin lives and a = x,y, z. These basic observables obey the 
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relations 




(1) 
(2) 
(3) 



For systems with a finite number of spins the energy of the system is 
encoded in the Hamiltonian and physically relevant states are ground and 
thermal states. The first are determined by the eigenstates of the Hamil- 
tonian corresponding to the lowest energy and the second are given by the 
Gibbs canonical density matrix at temperature T. The limit T — > of the 
Gibbs state is a ground state. In degenerate cases this limit is the average of 
all ground states and there are therefore generally more ground states than 
just the low temperature limit of the canonical Gibbs state. 

Typically one is interested in quantities that have a nice scaling behaviour 
with respect to N and hence in the leading asymptotic behaviour. This is 
precisely what the formalism of systems with an infinite number of particles 
catches: it provides right away a framework that allows to define and compute 
the asymptotic ground and equilibrium states. Obviously the total energy of 
an infinite systems is not a sensible observable and one characterises ground 
and equilibrium states in terms of the interactions between the spins. In 
fact sufficiently local interactions generate a well-behaved dynamics on the 
infinite system given by a group of automorphisms {a t \ t G R} of the 
observables. The generator of this group is a derivation 5 determined by 
a t = exp(it8). It satisfies Leibniz's rule S(XY) = S(X)Y + X 5(Y) and is 
self-adjoint S(X^) = — (5(X)) . For a finite system with Hamiltonian H one 



There are various ways of introducing ground and equilibrium states of an 
infinite system, such as thermodynamical limits of local Gibbs states possibly 
with boundary conditions or as states maximising the entropy density for a 
given internal energy density. We shall use the characterisation in terms of 
the dynamics. This leads to the ground state equation and the KMS condi- 
tion. For finite systems these conditions return the usual notions of ground 
state and canonical Gibbs state at a given temperature. 

A state oj satisfies the ground state equation if 



has S(-) = [H, •]. 




(4) 
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It satisfies the KMS condition at inverse temperature (3 = 1/ k^T > if there 
exists for any choice of observables X and Y a function 

F x ,y ■ {z e C | < 3m < (3} -> C (5) 

such that F X y is analytic in the open strip {z e C | < 3m < f3}, is 
bounded in the closed strip, extends continuously to the lines ^smz = and 
3m = (3 and satisfies for t G R 

F XtY (t) = oo{Xa t {Y)) and F xy {t + i(3) = oo{a t {Y)X). (6) 

Loosely speaking, the KMS condition says that 

uj(Xa^(Y)) = uj(YX) (7) 

at least if the expression a^Y) makes sense. It therefore links equilibrium 
states with the Heisenberg dynamics at imaginary times. For finite systems 
when the dynamics is given by 

a t (X) = e itH Xe~ itH (8) 

it follows readily from (7) that only the canonical Gibbs state at inverse 
temperature (3 satisfies the KMS condition. 

If the system has a local symmetry in the sense that there exists a local 
unitary U, i.e. a unitary living on a finite number of sites, such that 

at(UXrf) = Ua t {X)U ] (9) 

and if u is a KMS state, then also the transformed state X i— > uj(U X W) is 
a KMS state at the same inverse temperature (3 and is unitarily equivalent 
to the original one. A general property of KMS states, see[16], implies that 
uj is actually invariant under the symmetry 

lo(X) = u{UXU ] ). (10) 

There is no such result for ground states as determined by equation (4). 

We shall consider two models: the Ising ferromagnet and a model intro- 
duced by Kitaev. Both models have a finite version on a ring and a torus 
respectively. We shall also consider the infinite versions, they live on a chain 
and on a staggered square lattice. The particle interactions are in both cases 
nearest-neighbour and both models, either finite or infinite, allow for ground 
states in a very strong sense in that the states minimise all local interactions. 
This can be seen as a complete absence of frustration. 
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2.1 ID Ising ferromagnet 



The interaction for an Ising ferromagnet is given in terms of bonds 

Z b :=<r*tf +1 , with&={j,j + l}. (11) 

In this expression j is either a site of a ring or a chain. 

Let us first consider the finite system of N spins on a ring, meaning that 
we identify the sites 1 and N + l. Clearly the bond spins are not independent 
as they satisfy the relation 

n^ =i - ( i2 ) 

b 

Any (JV— 1) independent bond observables — one can omit an arbitrary one 
due to (12) — correspond to the a z observables of (N—l) qubits. The algebra 
A z generated by the bonds is Abelian but is not generated by a maximal set of 
commuting observables in the algebra of the N spins. In fact its commutant 
A' z within the spin algebra is the product of a qubit algebra and A z . This 
can be seen by the following direct construction. Choose 

X = o*o* ■■■<j x N , Y = a\al ■■■a x N , and Z = a{. (13) 

The matrices X, Y and Z satisfy the same relations as the Pauli matrices and 
commute with the algebra of the bond observables. It can be checked that 
any element commuting with A z belongs to the algebra generated by the 
observables of the encoded qubit (13) and A z . 

Let us now turn to the infinite system which lives on a chain. The local 
observables are finite linear combinations of finite products of single site 
spin observables and the quasi-local observables are obtained by completing 
this local algebra with respect to the natural norm. The completion of the 
algebra of the bond observables A z is again commutative but this time the 
commutant within the quasi-local algebra is also Abelian, it is generated 
by e.g. af and A z . The reason for this simplification is to be found in the 
disappearance of the infinite product of o x observables in the limit iV — > oo 
for the X and Y observables in (13). From a mathematical point of view 
such a limit does not exist in the spin algebra, from a physical point of view 
the distribution of such observables has no stable limit in any physically 
reasonable state. 
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We shall examine the ground and equilibrium states of the Ising Hamil- 
tonian with local Hamiltonian 

bcA 

Here A is any finite subset of the chain. 
2.1.1 Ground states 

To compute the ground states we can use the following property of a state: 
suppose that X — X< such that X < 1 and that oj{X) = 1, then uj{YX) = 
lo(XY) = uo(Y) for any Y. Indeed, < 1 — X and so we can extract the 
square root of 1 — X and apply the Cauchy-Schwartz inequality 

\uj{Y{1 - X))\ 2 < lo(Y(1 - X)Y ] ) oo(l - X) = 0. (15) 

Similarly co(XY) = u{Y). 

Suppose now that w is a state on the Abelian algebra A z generated by 
the bond observables Zf, such that uoiZh) = 1 and let Co be an extension of 
this state to the spin algebra. As Zb < 1, we can apply the argument of 
above and conclude that for any element X in the spin algebra and any bond 
b£o{XZ b ) =u{Z b X) =u{X). 

This implies that Co is a ground state in the sense of (4). It suffices to 
show that Co(X*5(X)) > for any local observable X. For such observables 
the action of the derivation 5 is given by a finite sum of commutators 

S(X) = -Y^[Z b ,X\. (16) 

6eA 

Here A is the set of bonds which have a non-empty intersection with the 
dependency set of X. We now compute 

Co(X ] 5(X)) = ^Co{X ] X Z b - X ] Z b X) (17) 

6eA 

= ^(^(X f X) -Co(X ] Z b X)^j (18) 

beA 

= ^2<j(X*(l-Z b )X) > 0. (19) 

6eA 
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To compute the ground state expectation of an arbitrary spin observ- 
able for the finite system we consider arbitrary monomials in Pauli matrices. 
Using 

u(o?) = -u{Z b <j* Z b ) = -Cu(aJ) (20) 

for a site j which belongs to the bond b and the obvious extensions to arbi- 
trary polynomials, we conclude that only monomials belonging to the corn- 
mutant of A z can have non-zero expectation values. Moreover with X as 
in (13) 

Cu(XB) = u(X)cu (B), B e A z , (21) 

and similar expressions replacing X by Y or Z. The state ujq is the unique 
state on A z that assigns the value 1 to all products of bond observables. 
We have therefore determined all ground states of the finite model: they are 
products of an arbitrary state on the qubit algebra and the state u . 

The infinite system can be handled in exactly the same way, we now 
obtain product states of a classical bit with the state u . In particular, this 
implies that the set of ground states of the Ising ferromagnet on an infinite 
chain is a simplex, consisting only of mixtures of two extremal states. In 
principle, this could be used to encode a classical bit. In fact, this is a bit 
too optimistic as it is well-known for this model that at finite temperatures 
there is a unique equilibrium state. 

2.1.2 Equilibrium states 

Determining the equilibrium states of the finite and infinite Ising models as 
solutions of the KMS condition can be carried out along the same lines as in 
the previous section. Using for any site j belonging to a bond b 

Z b a*Z b = -a* and Z b o)Z b = -o) (22) 

and the general result about KMS states mentioned in (10) we conclude 
that the only monomials that have non-zero expectation values belong to the 
commutant of A z . The solution of the KMS condition is the extension of the 
equilibrium state for the bond observables to A' z assigning zero values to any 
observable of the form XB, YB or ZB, B e A z for the finite system and zero 
value to o\B for the infinite system. 
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2.2 Kitaev's model in 2D 



The finite version of Kitaev's model lives on a K x K lattice with periodic 
boundary conditions (toroidal lattice), see the square dots in Fig. 1. The 
spins live on the edges of the lattice as shown by the black and white dots 
in the figure and so we have N = 2K 2 qubits. The interactions between 
the spins are given by star and plaquette terms X s and Z p . A star is a 
diamond whose vertical vertices lie on the lattice of black dots, while the 
vertical vertices of a plaquette lie on the lattice of white dots. The actual 
observables are then 

X s = Y[a? and Z p = \\a) (23) 

jes jep 

The names star and plaquette correspond to the squares and stars marked 
by thick lines in the figure. 



HI O- 

II 

HI O- 

II 

HI O- 

II 

Fig 1: Kitaev's lattice 

As stars and plaquettes have either or 2 sites in common, [X s , Z p ] = 0. 
So, the algebras A x and A z generated by the X s and Z p are Abelian. We shall 
denote by A xz the algebra generated by A x and A z . Because of the periodic 
boundary conditions the star and plaquette observables are not independent, 
they satisfy the relations 

J]X S = 1 and n^ = 1 ' ( 24 ) 

s p 

10 




otherwise they are independent. As in the Ising case the algebra A xz is 
not maximally Abelian within the spin algebra. Its commutant consists of 
a product of two qubit algebras and A xz . This can again be seen quite 
explicitly by the following direct construction. Similarly to (13), one can 
introduce observables for two encoded qubits 



jeci jec 2 
Here ci, c[, ci and c' 2 are the loops shown in Fig. 2. 



(25) 



o 



o 



-©- 



f2 



Fig 2: The loops in the toric lattice 



Let us now turn to the infinite system which lives on an infinite staggered 
lattice. In contrast to the Ising case the commutant of the algebra A xz in 
the spin algebra is now A xz itself. This can again be understood in terms of 
the non-existence of the limits of the loop operators in (25) when iV — > oo. 
Indeed with increasing N, the length K of the loops increases as N 1 ^ 2 . 

The local Hamiltonians of the model are 

#A itacv = - E X °- E Z v ( 26 ) 

{s|sCA} {p|pCA} 

Here A is any finite subset of the staggered lattice. 
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2.2.1 Ground states 



Let u be the state on the Abelian algebra A xz generated by the A x and A z 
such that 

u (X s )=u (Z p ) = 1. (27) 

Using the same arguments as for the Ising model we obtain that ujq ex- 
tends to a product state of a two qubit algebra and A xz for the finite model. 
For the infinite model there is only a unique extension Q as we argue in the 
following lines and this has, by unicity, to be pure. As before, we compute 
the expectation values of a finite product of elementary spin observables us- 
ing (15). Consider a product of Pauli matrices living at the sites of a finite 
set A. Consider first the point of A farthest to the North and also, on the row 
to which it belongs, farthest to the East. Suppose that this point belongs 
to the lattice of black dots. An argument similar to that in (20) kills the 
expectation of the observable unless it carries a o z Pauli matrix. Suppose 
it does and consider then the point on the lattice of dots, lying South-East 
of the first. Unless this carries also a a z spin the expectation also vanishes 
because of anticommutation with a suitably chosen plaquette observable. If 
both points carry such a z matrices, then they can be removed by multiplying 
with a plaquette observable having the two points as North-East boundary. 
A similar argument holds if our original point belonged to the lattice of white 
dots. Continuing like this we see that the only observables having non-zero 
expectation value are products of star and plaquette observables, but all 
these have expectation 1 because of (15). So this model has a unique ground 
state and cannot be used to encode even a single classical bit. 

It remains to show that indeed a state satisfying (27) exists. It suffices 
therefore to observe that we can map the star observables on spins of an Ising 
model on Z 2 and the plaquettes on these of a similar Ising model. Doing so, 
the interaction in the Kitaev model just becomes that of a couple of free 
Ising spin systems on Z 2 , which clearly admits such a ground state. 

2.2.2 Equilibrium states 

Probably the most important feature of the plaquette and star observables 
is that 

[X S ,X S /] = [X s , Z p ] = [Z p , Zp>] = (28) 
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for all choices of s, s', p and p' . Moreover, 

X a = Xl Z p = Zl and X s 2 = Z 2 p = 1. (29) 

Therefore, for the infinite system, the algebra generated by the X s and Z p 
is isomorphic to the algebra of continuous function on a configuration space 
which is a countable product of copies of Ising spaces {til}- We can, more- 
over, choose the isomorphism in such a way that each of the X s and Z p is 
mapped onto an Ising variable. 

Let A C Z 2 be a finite set. Adding successively boundary layers, we 
obtain a collection of sets 

A(0) := A (30) 

A « : =( U s )U( U p) ( 31 ) 

s nA(o)^0 P nA(o)^0 

A(2):=( [J *)U( U p) (32) 
s nA(i)^0 p nA(i)^0 



In the Heisenberg picture, a local observable X initially living in a finite set 
A evolves after time t into 

a t (X) :=X + itS(X) + ^S 2 (X) + ■■■ (33) 

where the derivation 5 is obtained from the local Hamiltonians 

6(X)=lim[H™~",X] (34) 

M— >Z^ 

= - [*.>*]- E (35) 

sCA(l) PCA(1) 

= [Ef^\X\. (36) 
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Using the commutation relations between the X s and Z p , we see that 

= ~ E - E ( 37 ) 

scA(2) pCA(2) 

= " E E [*.,[*✓,*]]"••• (38) 

scA(2) s'cA(l) 

= " E E [X.s',[Xs,X]}-... (39) 

s'cA(l) scA(2) 

= " E E lXs>,[Xs,X}]-... (40) 

s'cA(l) sCA(l) 

We therefore have for any A: = 1,2,... 

5 k {X) = [H^f v , [Hf$ ev , ■ ■ ■ [H^f v , X] ■ ■ ■ } } k commutators. (41) 

Hence, for all i G R, a t (X) lives in A(l). 

The commutation of the X s and Z p also allows us to conclude that for 
any s and p and any quasi-local observable X 

at([X a ,X]) = [X s ,a t (X)] and a t ( [Z p , X]) = [Z p , a t (X)] . (42) 

Because of (42), given an a t -KMS state u, for any choice of s and p and 
for any s£R, also 

X ^ iu(e isXs Xe- isX ^ and XM W ^ sZp le^) (43) 

are a^-KMS states at the same inverse temperature f3. All such states are 
unitarily equivalent to the original one. By (10), uj is invariant under all 
these local unitary transformations 

lu(X) = uo (e isXs Xe~ isxA = oo (e isZp Xe~ isZ ^ . (44) 

We now compute the ctj-KMS states. It suffices to compute the value of 
uj on any monomial in the Pauli matrices, i.e. on any observable of the kind 

X = W (T 7 e,e{x,y,z}. (45) 
jeA 
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Using the invariances (44) with the value s = n, we find that 

u>(X)=u,(j[v?)=0, (46) 
jeA 

unless X is a product of star and plaquette observables. This implies that uj 
is completely determined by its values on the commutative algebra generated 
by the A x and A z . The restriction of uj to this algebra is, however, a classical 
equilibrium state of a free Ising Hamiltonian at inverse temperature f3. It is 
well-known that there is only one such state which has a product structure 
fully determined by the expectations 

uj(X s ) = uj(Z p ) =tanh(|). (47) 



3 The open Ising and Kitaev models 

We shall for both models consider noise originating from a coupling of the 
atomic spins to a surrounding heat bath. In order to keep things as simple 
as possible, we shall assume that each spin has its private heat bath which is 
moreover independent from these of the other spins. All baths are assumed 
to be identical. Moreover, we shall restrict ourselves to the Markovian ap- 
proximation applying directly the weak coupling limit of Davies [6]. 

Before treating the Ising and Kitaev models, we briefly sketch the general 
setup and properties of Davies generators. A small system is coupled to a 
collection of heat baths leading to the global Hamiltonian 

H = H sys + H bath + H' mt with H' mt = J2 S «® f«, (48) 

a 

where the S a are system operators and the f a bath operators. The main 
ingredients are the Fourier transforms h a of the autocorrelation functions of 
the f a . The function h a describes the rate at which the coupling is able to 
transfer an energy difference uj from the bath to the system. Often a minimal 
coupling to the bath is chosen, minimal in the sense that the interaction part 
of the Hamiltonian is as simple as possible but still addresses all energy levels 
of the system Hamiltonian in order to produce finally an ergodic reduced 
dynamics. The necessary and sufficient condition for ergodicity is [26, 27] 

{S a , H^}' = Ct, (49) 
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i.e. no system operator apart from the multiples of the identity commutes 
with all the S a and # sys . 

We begin by introducing the Fourier decompositions of the SaS as they 
evolve in time under the system evolution 

e itH°y° Sa e -itH°y° = J2 Sa (u) e~**. (50) 

U! 

Here the to are the Bohr frequencies of the system Hamiltonian. From self- 
adjointness we have the relation 

S a (-u) = S a (u)l (51) 

The weak coupling limit procedure then returns the following equation for 
the evolution of the spin system in Heisenberg picture 

^ = ziH^, X] + C dis (X) =: C(X) (52) 

For thermal baths one has moreover the relation 

h a (-u)=e- pu h a (u) (54) 

which is a consequence of the KMS condition. The operator C generates a 
semigroup of completely positive identity preserving transformations of the 
spin system. It describes the reduced dynamics in the Markovian approxi- 
mation and enjoys the following properties 

• The canonical Gibbs state with density matrix 

e -/3H s y s 

Pp = 7 V ( 55 ) 

Tr(e-^ sys J 

is a stationary state for the semigroup, i.e. 

Tr( P/3 e t£ (X)) =Tr{p p X). (56) 
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• The semigroup is relaxing, meaning that for any initial state p of the 
system 

lim Tr(pe tc (X)) = TrippX). (57) 

• Furthermore, the generator satisfies the detailed balance condition, of- 
ten called reversibility Writing S(X) := [# sys ,X], 

[5, £ dis ] = and Tr ( P/3 yt £ dis (X)) = Tr( P/3 (£ dis (y)) f x) . (58) 

The last equation expresses the self-adjointness of the generator with 
respect to the scalar product defined by the equilibrium state. The 
space of observables equipped with the scalar product 

(X,y) /3 :=Trp /3 Xty (59) 

is called the Liouville space and the generator of the reduced dynamics 
is a normal matrix on that space, i.e. the Hermitian and skew-Hermitian 
parts of the generator commute. 

3.1 Dissipative generators for Ising and Kitaev models 

The eigenstates and eigenvalues of the Ising Hamiltonian are conveniently 
labelled by the excited bonds and one more dichotomic variable which takes 
the twofold degeneracy into account. This last could be the eigenvalues of 
Z = o\. An excited bond is often called kink and can be seen as a quasi- 
particle. In the model on a ring, kinks always appear in pairs due to the 
topology of the system. The two dimensional Kitaev model will admit a very 
similar description in terms of excited plaquette and star observables now 
called anyons. 

More precisely, for the Ising model on a ring, the eigenstates of the Hamil- 
tonian are written as |&i, b 2 , . . . , b 2n ; z), they have energy 

# Ising |6i, 6 2 , . . . , b 2n , z) = (2n - N) \h, 62, ... , b 2n , z) (60) 

The interaction Hamiltonian is chosen as 

N 

ir" ;>>•; ./• (6i) 
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where fj is a self-adjoint field of the j-th bath and all the baths are isomor- 
phic. 

The Bohr frequencies of the Ising Hamiltonian are 0, ±2, ±4, .... Only 
the Bohr frequencies and ±2 contribute to the Davies generator (53) due 
to the choice of the interaction. We compute the Fourier components of the 



ofs 



e i«PP* a x e -itH^ = Q -2it a . + Q 2it Q f + Q (62) 

The operators dj and a° are given in terms of projection operators 

P? = \ (1 - Z h Z v ) and P± = 1 (1 =f ^) (1 =F ^) (63) 
with j = bn b'\ 

«; p >* p ? p >* ^ /, ; , ( 64 ) 

%- / 7^' / 7- / 7^-^' / 7- ( 65 ) 

The Davies operators allow for an interpretation in terms of kinks. Consider 
two bonds b and b' having a site j in common. The operator a° kills states 
for which both b and b' are either empty or occupied and exchanges the 
occupations if only one of them is. The operators and a~ create and 
annihilate states with two kinks meeting at j respectively and kill the other 
states. This leads us to the master equation in Heisenberg picture 

N 



f \ E{^( 2 ) (°5 ^ °j] + % + e ~ 2/3 % [*> °J] + e ~ 2/3 K' x ] a 
-£(o) [ajjaj,*]]} (66) 



For Kitaev's model on a torus, the eigenstates of the Hamiltonian are 
written as \pi,P2, ■ ■ ■ ,P2m', «i, S2, • • • , s 2n ; zi, z 2 ). The labels z\ and ^ 2 refer to 
the eigenvalues of the z-components of the encoded qubits while the pj and 
s k point to the excited plaquette and star observables. Such an eigenstate 
has energy 

H Kltacv \pi,p 2 , ■ ■ ■ ,P2m] si, s 2 , • • • , s 2n ; zi, z 2 ) 

= (2m + 2n- N) \p u p 2 , . . . ,p 2m ; s u s 2 , . . . , s 2n ; z u z 2 ) (67) 
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The interaction Hamiltonian is chosen as 

N 2N 
3=1 j=N+l 

where fj is a self- adjoint field of the j-th bath and all the baths are isomor- 
phic. 

The Bohr frequencies of the Kitaev Hamiltonian are 0, ±2, ±4, . . .. 
Again, only the Bohr frequencies and ±2 contribute to the Davies gen- 
erator (53) due to the choice of the interaction and we compute the Fourier 
decompositions of the spins as they evolve in time 

e itH^ a x e -itH^ = e -2it a . + e 2it fl t + fl (6Q) 

e itJJ8pi ° o) e- UHSpia = e~ 2U b 3 + e 2U b) + 6°. (70) 
The operators a and b are again given in terms of projection operators 

P* = \{t-Z p Z p ,) and Pf= l -{ 1T Z P ) (1TZ P ,) (71) 

with j = p(~)p'. We need a second set of projections associated with the stars 

= l - (1 - X S X S ,) and i?± = l(l T X s ) (l =f (72) 

with j = s n s'. Then 

^-'f^'f-'f^-^'f (73) 

/' < /y /' < /y • (74) 

and 

b° = R°a*R° = R°a* = a*R° (75) 

^ / ' ; °j R j (76) 

The Davies operators allow for a similar interpretation as in the Ising case 
in terms of anyons. This leads us finally to the master equation in Heisenberg 
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picture 

r rrKitacv v i 
N 

1 , 

+ 



^{^(2) (at [X, aj ] + [at, X] a 3 + e" 2/3 a, [X, at] + e~ 2/3 [a„ X] a}) 
-MO) K°,[a°,X]]} 

fc(O) (77) 
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l , 

+ 



Both the Davies generators for the Ising and Kitaev models extend to 
generators of completely positive identity preserving semigroups on the in- 
finite chain and the infinite staggered lattice. This easily follows from the 
locality of the elements in the generator and the non-propagation of both the 
Hamiltonian and dissipative parts in the generators. 



3.2 Dynamics of the qubit observables 

In order to treat both the Ising and Kitaev models in a unified way it is 
useful to introduce the notation Aab for the Abelian algebras of observables. 
In the Ising case A^ = A z and for the Kitaev model ^4 ab = A xz . In both 
cases we denote by C the commutant of ^4 ab in the respective spin algebras, 
so C = QA Z for Ising and C = Q\ <2 2 A xz for Kitaev. The algebras Q and 
Qi 5 2 are the qubit algebras introduced in (13) and (25). So, e.g. in the Ising 
case any element of C can be written as linear combination of 1 and the qubit 
observables X, Y and Z with coefficients in A a b- Let us denote generically by 
{Q} the bases of the one or two-qubit spaces, so {Q} = {l,X, Y, Z} in the 
Ising case. 

We consider first more general qubit observables {Q} = {1, X, Y, Z} which 
enjoy the following stability property: [Q, A] = for all A e *4. a t>, i- e - {Q} c 
C. The physical meaning of this condition is that the {Q} are not only 
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invariant under the Ising or Kitaev dynamics but also under a whole class of 
perturbed dynamics of the form 

Rising ;= _ J b Z b and ij'Kitacv . = _ ^ jx ^ _ ^ jz ^ ^g) 

In these formulas, the coupling constants are arbitrary but strictly larger than 
0. Remark that these perturbed models return precisely the same ground 
states as the original ones. In fact it will suffice to consider encodings of the 
type 

Q = QF Q , F q e (79) 

where {Q} is the choice mentioned above. In order to satisfy the qubit 
relations we have to impose on the functions F that 

F X F Y = F Z , F Q = Fl, F% = 1, and Ft = 1. (80) 

In order to quantify the quality of a memory we are looking for the short- 
est decay time of the time autocorrelation functions 

(Q(f),Q>/3 = Tr( P/3 e t£ (Q) q) = Tr( P/3 e«^ (Q) q) = \\Q(t/2)\\% (81) 

This chain of equalities follows from the properties of the generator of the 
reduced dynamics, derived at the beginning of Section 3. The time auto- 
correlation functions measure how well a system encoded with a message at 
time has retained its information up to time t. The rightmost Q in the 
product can be seen as a polarisation of the system corresponding to writing 
a message at t — while the leftmost is the readout at time t. 

We shall now compute the action of the reduced dynamics on elements 
of the form Q = Q Fq with Fq as above and show that 

e tc (Q) = Q e'^ (F Q ) = e*^ (F Q ) Q (82) 

with 

e* £ ~ Q : Ab -> Ab (83) 
a contracting semigroup, hence 

(Q(t),Q) /3 = <e^(F Q ),F Q > (84) 
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In order to verify (82) we need some relations 

a? F E A ah whenever F e A ah (85) 
o?Qo? = g%(j) Q with e {±1}- (86) 

Consider a typical contribution to the generator 

K(X):= aj X a] - \a ja ) X - \X aj a) (87) 

with 

aj = Pra?P+, P±eA ah . (88) 

For F G Aah we obtain 

K(Q F) = PJ o] Q o] cr| F <rj P~ -QP~ F (89) 
Q(/ , ; (.'/Q(./)rr;/-n : ; F)) =:QKq(F). (90) 

Clearly £ Q (F) e A ah . It now suffices to add all contributions in order to 
show (82). In this way we obtain the following expressions for the reduced 
generators Cq 

N 



Rising (F) = J- D» (<£(j) (Tj F<Jj — F) (91) 

J'=l 
AT 

Z K it acv (F) = £ ^) aJFaJ - f) 
j'=i 

+ ^(sqOK^-^)} (92) 

F>* = fe(-2) P- + fc(2) P+ + h(0) P° (93) 
£>? = h{-2) Rj + h(2) R+ + h(0) P°. (94) 

Finally, let us analyse the evolution of some special observables in Kitaev's 
model. Consider observables Xi and Zi of the first qubit given by 

Xx = X ± F X , Zi = Z 1 F Z , (95) 
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with F x G A x , F z G A z , and Xi, 7-\ as in (25). It suffices to consider one of 
them, e.g. Z 1 . It turns out that the following Hamiltonian and coupling with 
the environment generate the same time evolution of this observable 

TV 

^Kitaev = _ Z p and H' mt = J2tf® fj- (96) 
p j=i 

More precisely, we have 

*Z i = e £— *z i (97) 

where £ Kltacv ^ i s solely determined by Davies operators of a-type 

<<] °1 P * and «i ~ °V7- ( 98 ) 

Similar relations hold for the second observable. This can be interpreted as 
follows: with this choice of basic observables, bit and phase evolve separately 
according to the same evolution up to an exchange of o x and a z . This 
property is then inherited by the reduced generators. Such a separation is 
not surprising. It simply means that such a choice of qubit observables is 
compatible with the structure of the quantum code behind Kitaev's model. 
Namely, the subspace of ground states constitutes a so-called CSS code [28, 
29]. The characteristic feature of this class of codes is that the recovery 
procedure, aiming at restoring the initial state of the qubit after the attack 
by the noise, is divided into two stages: first correcting bit errors and then 
phase errors. The correction procedures are again the same, modulo a z <-> a x 
exchange. 

The aim of this paper was to analyse the static and dynamical structure 
of Kitaev's model in 2D, using the analogy with the Ising ring or chain. 
More involved models in 3 and 4D should be addressable in a similar way. A 
further topic is determining the spectral properties of the reduced semigroup 
generators, especially the behaviour of relaxation times as a function of the 
system size. Such questions are important in order to decide whether such 
systems are reasonable candidates for stable quantum memories. This will 
be dealt with in a forthcoming paper. 
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